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Abstract. Waves associated to large scale oceanic motions are gravity 
waves (Poincare waves which disperse fast) and quasigeostrophic waves 
(Rossby waves). In this Note, we show by semiclassical arguments, that 
Rossby waves can be trapped and we characterize the corresponding 
initial conditions. 



1. Introduction and results 

Large scale oceanic motions can be described - at first sight - by the 
linearized Saint- Venant equations for thin layers, with Coriolis force, i.e. by 
the following two-dimensional system [6]: 

(1) 

(0 d xi d X2 \ / p 

d Xl -B(x 2 ) J7 = 0, where U = U(t, x u x 2 ) = ui 

8 X2 B(x 2 ) J \u 2 

where B is the local vertical component of the Earth rotation vector (de- 
pending only on the latitude x 2 ), p denotes the fluctuation of water height 
and u the horizontal velocity field. For the sake of simplicity, we assume 
that (xi,x 2 ) G M x T, meaning that we neglect the influence of the lateral 
boundaries. 

Waves associated to that linear system are usually classified in two 
families : First, Poincare waves, which are fast dispersive gravity waves; 
Second, Rossby waves, due to the inhomogeneities of B, which propagate 
much slower [2], [3]. 

In this Note, we are interested in the propagation of Rossby waves. We 
show, in agreement with physical observations [TJ, that these waves can be 
trapped in some regions, called ventilation zones, which are not influenced by 
external dynamics and sources such as continental recirculation for instance. 

More precisely, in the limit of large values of B = b/e, we construct 
a co-dimension 1 submanifold A of T*(M x T). This set contains the e- 
microlocalization region inside which, essentially, an initial condition of (pQ) 
remains trapped. Let us first recall that the e- frequency set of a function u 
[I] is the (closed) subset of T*(M x T), complement of the set of points (x, £) 
such that there exists a C°° function x, with = 1> such that 
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Theorem 1.1. Let us consider the system 

U = with an initial condition 
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d xi 


dx 2 


(2) d t U + d xi 





b(x 2 ) 

€ 


\ d X2 


6(3:2) 

e 






(3) C^| t =o = «5l 



mi/i L 2 (M x T, dx\dx 2 ) conditions (periodic inx 2 ). Let A = {F(£\, x 2 , £2) 
} c T*M x T, where F is defined in Lemma \3.1\ 

We suppose that the e-frequency set of U\t=o is contained in a compact 
set C satisfying: 

(4) c n {6 = 0} = c n + Kx 2 f = 0} = 0. 

Let us fix a compact set Q o/M x T. Then there exists (e-)pseudo-differential 
operators P9, P®, P® of principal symbols p Q p = ib{x 2 )^i{^2 + £1 + b 2 {x 2 ))~ l , 
P°i = -66(£l + + b\x 2 ))-\ p\ = eM + & + b\x 2 ))-\ such that: 

(1) if the e-frequency set of PpP° + P^u\ + P 2 u 2 does not intersect A n 
T*n, then 3T > such that: 

\\U(- e )\\ LHQ) = 0(e^)fort>T. 

(2) if the e-frequency set of P p p° + P\u1 + P 2 u 2 does intersect A n T*f2, 
then, \/t > 0, \\U(j)\\ L 2(ty ^ 0(e°°) (in other words the frequency 
set ofU(-) intersects T*Q). 

In the case of a WKB initial condition the conclusion is more precise. 
Consider 

R\x) s _ 
U\ t=0 = V?(x) I e*^ 



X 



Suppose than the (Lagrangian) manifold 



{(x,VS(x)), (p° p R° +p\U% +p 2 ) U$)(x,VS(x)) ^0} intersects AnT*ft, then: 



(5) \\U(- e )\\ L2{n) ~C(t) + 0(e), C(t)>0. 

2. Reduction to a scalar situation 

Performing first a Fourier analysis in x\ (as the system does not contain 
explicitly this variable) and looking secondly at modes of the system ([2]) , we 
can prove the following Proposition, heart of our results: 

Proposition 2.1. There exist three pseudo- differential operators T±, To, of 
leading symbols 
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r± = ±T/£i + & + b 2 (x 2 ), T = e $+$+jp\ X2) , roots of fW, such that, if 
is microlocalized outside £f — rf = 0, we have, for j G {— , 0, +}: 

{ied X2 T ] +ed Xl b{x2)){e 2 dl 1 +Tj)- 1 
ed t u 2j = iT jU2iJ => Uj := ( -(e 2 8 Xl d X2 + ib(x 2 )T j )(e 2 d 2 Xl + if)- 1 ) u 2) 

Id 

satisfies up to 0(e°°). 

The proof is a consequence of the result contained in the Appendix. 
The following result shows that any initial condition of © can be decom- 
posed on the three modes of the last section. 

Proposition 2.2. \/p,u\,u 2 , 3u 2 j, j G {+,0,—} such that: 
(6) 

p\ ( (ied^Tj + ed^bfam^+Tf)- 1 

111 = E -(^ X1 d m + ibix^T^dl + If)- 1 | u 2>j + 0(e c 
u 2 j Id 



■-'•J2Q?U2j + 0{e°°) 



j 



To prove the Proposition one has just to invert the matrix: 
Semiclassically it is enough to show that the matrix of the leading order 
symbol is invertible 



(7) 



£li2+ib(x2)^/( 2 +i 2 '+b(x 2 ) 2 _^ Slg2-»fc(^2) V /Sf+ej + t'(3:2) 2 

e 2 +Hx 2 y 2 ~«T if+f^y 2 

\ i i i j 



A simple computation shows that the jacobian J = ^gt^^J^^ — — 2- In 

particular the inversion of the matrix ((QT Q° Q+) can be done symbolically 
at any order and gives the leading order symbols of the operators P J := 
(F%,P[,P$) such that u 2j = P 3 p p + P{m + P 3 2 u 2) j G {-,0,+}. One gets : 
P° P = ib{x 2 )tM + & + b 2 {x 2 ))~\ p§ = + + b 2 ^))- 1 , p§ = 

eM+ei+b 2 ^))- 1 . 

3. Propagation under the Rossby Hamiltonian r 

Thanks to Proposition 12.11 it is enough, as far as the Rossby mode is 
concerned, to look at propagation with respect to the Hamiltonian tq/e. 
This Hamiltonian being independent of xi, £i will be conserved. The flow is 
periodic in the variables x 2 , £ 2 (one degree of freedom). Therefore, since x\ = 

b ^r^rp^7~Tp^ is periodic (the case of infinite and zero period is treated 

in [1]), x\(t) will contain a part, linear in time except if 



f(p r rm f tow - r ri ° d ^))M)!^l±j!(^M ^ - n 
*■(&,**«>), & (o» - y o (6(^T^R^ * " °' 
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Lemma 3.1. As b'(x 2 ) + 0, V {x 2 )F{^ x 2 , £ 2 ) >0os(i 
£i ^ and is invariant under the flow of ^ . 

Define Efa,x 2 ,Z 2 )- b ' {x ^ 



±oo, < as 



Then |-F(£i, zr 2 , £2) 



x+ 



b'(x) 



2£i 



E(6,x a ,&) 



i dx 



et-b^x) 



where 



zs i/ie largest interval of T containing x 2 in which Wg^^g^ 

We define A := {(x ls a^fr, 6)/ F(^,x 2 ,^ 2 ) = 0}. 
Thanks to Lemma 13.11 A ^ and dimA = 3. 



6 2 (x) > 0. 



Corollary 3.2. Suppose b'(x 2 ) / 0. Then \F(£x,x 2 ,£ 2 )\ > as ^ 

0, with C(x 2 ,^ 2 ) > 0. This implies that the trapping phenomenon will take 
place only with initial conditions oscillating enough in x\. 



Remark 1. : Since the Hamiltonian E does not depend on x\, we can ex- 
press it, for each value of £i, on the action variable A: E(£i,x 2 ,£ 2 ) = 
H(A,£i). This allows to define the function A(£i,x 2 ,£ 2 ) by E(£i,x 2 ,£ 2 ) = 
H(A(£i,x 2 ,£ 2 ),£i). One can easily show that: 

F(6 ' X2 ' 6) " d A H(A,^) U = A ^ x ^ 

and the following variational characterization of A: Let us fix the energy to 
E and let T(^i,E) be the energy shell {E(^i,x 2 ,^ 2 ) = E}. Then 

A={j{(m,E),&}, 

E,i 

where £'[ is such that the area inside T(£\,E) is extremal. 

4. Dispersion of Poincare waves and proof of the Theorem 

The proof of the Theorem involves, for the Rossby modes (j = 0), the 
standard result of propagation of the frequency set. If the initial frequency 
set is such that part of the Rossby mode is trapped, in particular if it 
intersects {b'(x 2 ) = 0}, this concludes the proof. 

If not we have to prove some dispersion for the Poincare modes (j = ±) for 
times of the order = , for which the theorem of propagation of the frequency 
set is not enough. But, the system being integrable, we can perform an 
expansion on (Bohr-Sommerfeld) eigenvalues of the Hamiltonians T± and a 
decomposition of U\t=o on a compact set of coherent states (thanks to the 
condition on the microlocalization of the initial datum) [5]. 

Since t± = ±vCl + b 2 (x 2 ) + £f + O(e) and we are microlocalized far 
away from £f + b 2 (x 2 ) + Ci = 0, we can find pseudo-differential operators 
H 2 ± of principal symbols £| + b 2 (x 2 ) such that T± = ±^H 2 ± + £f . The 
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Bohr-Sommerfeld quantization condition (with subsymbol) gives that the 
eigenvalues of Hz± are of the form: 



A± (W^)+e4(6) + 0( e 2 



where A± is the energy £| + b 2 {x2) defined on action variable, and e/z±(£i) S 
C°° is the correction due to the subsymbol. Propagating at time t = s/e a 
function, product of a coherent state at (q, p) (in xi) and an eigenfunction 
of T± (in X2), gives rise to expressions of the type: 

, e(xi - g)ft ± (A| + jf)W ze(fr - p) 2 

exp 1 a^i. 

e 

The stationary phase lemma then gives that this integral is 0(e°°) except if 
there exists a stationary point, given by the conditions: 

(2ft+e%/4) 
£1 = p and e(xi — gj ± — s = U. 

The second condition gives: 2yJ \± + £^ (xi — = -F + ^£iM±) s - There- 
fore, since p / and the X k 's are bounded by the above condition ([5]) on 
C, there is no critical point for x\ in a compact set. 

Appendix A. A microlocal analysis Lemma 

In this Appendix we give the crucial lemma for the proof of Proposition 
12.11 It tells us that the principal symbols of T±, Tq can be computed by 
solving the symbolic equation associated to 

2o2 ,2/ \ ■ 2^ { X 2)9 X ± 2 , 2o2 \ n 

e o; 2 - b (X2) - it — - + r + e J u 2 = 
obtained by (exact) algebraic computations from ([I]), that is 

(8) T 2_£2_£2_ 6 2 + 3l =Q 

r 

Lemma A.l. Let h = h(x, £, a) be a smooth function such that d a h\h = Q 7^ 0, 
and let r = t(x,£) be any continuous root of h(x, £, r(x, £)) = 0. 

Then there exists a pseudo-differential operator T = T(x, —ied x ) with 
principal symbol r(x, £) such that: 

(9) Tip = Xip H{X)iP = 0(e°°) 

where H(\) is a pseudo-differential operator of full symbol h(x,£, A). 

Sketch of the proof: the proof uses essentially pseudo-differential func- 
tional calculus [4]. 

For any (x, £) the principal symbol of h(x, £, T) at (y, £') is h(x, £, r(y, £')). 
We can write at first order: Hip := J e % ^ M ~h(x,^,T(y,—ied y ))ip(y)^ 1 = 
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J e i i -^ e i^ S ^h{x,i,T{y\i'))^{y')^^L. The integral over y gives <5(£- 

£')■ Therefore: Hij)(x) = /e^/i^.^rfi/.O^fe)^. So the principal 
symbol of H is h(x, £, t(x, £)) which, by assumption, is 0. 

For the e°° result, it is enough to repeat the same argument taking into 
account lower order terms e k Tk for and adding to any symbol a term of the 
form X^>i efc -f , fc T where the P^s are differential operators. We end up with 
an equation for r e ~ r + Yl ^ T k of the form: 

h(x, £, r e ) + J2k>i eh Qk{T, ■ ■ ■ , d l x d™T e ) = 0, that can be solved recursively 
under the condition d T h(x, £, T)\ h ^ x ^^ =0 ^ 0. 

Acknowledgements 

We would like to thank Agis Athanassoulis and Bach-Lien Hua for helpful 
discussions. 

References 

[1] C. Cheverry, I. Gallagher, T. Paul and L. Saint-Raymond, in preparation. 

[2] A. Dutrifoy, A. J. Majda and S. Schochet, A Simple Justification of the Singular Limit 

for Equatorial Shallow- Water Dynamics, in Communications on Pure and Applied 

Math. LXI (2008) 0002-0012. 
[3] I. Gallagher and L. Saint-Raymond, Mathematical study of the betaplane model: 

equatorial waves and convergence results. Mem. Soc. Math. Fr. (N.S.). 107 (2006), 

v+116 pp. 

[4] A. Martinez, An introduction to semiclassical and microlocal analysis, Springer (2002) 
[5] T. Paul, Echelles de temps pour revolution quantique a petite constante de Planck, 

Seminaire X-EDP 2007-2008, Ecole polytechnique. 
[6] J. Pedlosky, Geophysical fluid dynamics, Springer (1979). 
[7] J. Pedlosky, Ocean Circulation Theory, Springer (1996). 
E-mail address: christophe.cheverry@univ-rennesl.fr 

E-mail address: gallagher@math.jussieu.fr 

E-mail address: thierry.paul@ens.fr 

E-mail address: Laure.Saint-Raymond@ens.fr 

(Cheverry) Institut Mathematique de Rennes, Campus de Beaulieu, 263 av- 
enue du General Leclerc CS 74205 35042 Rennes Cedex 

(Gallagher) Institut de Mathematiques de Jussieu et Universite Paris VII, 
Case 7012, 2 place Jussieu 75251 Paris Cedex 05 

(Paul) CNRS and DMA Ecole Normale Superieure, 45 rue d'Ulm, 75230 Paris 
Cedex 05 

(Saint-Raymond) Universite Paris VI and DMA Ecole Normale Superieure, 
45 rue d'Ulm, 75230 Paris Cedex 05 



